Electrical control of magnon propagation in multiferroic BiFeOa films 
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The spin wave spectra of multiferroic BiFeOa films is calculated using a phenomenological Landau 
theory that includes magnetostatic effects. The lowest frequency magnon dispersion is shown to be 
quite sensitive to the angle between spin wave propagation vector and the Neel moment. Since elec- 
trical switching of the Neel moment has recently been demonstrated in this material, the sensitivity 
of the magnon dispersion permits direct electrical switching of spin wave propagation. This effect 
can be used to construct spin wave logical gates without current pulses, potentially allowing reduced 
power dissipation per logical operation. 



One of the challenges of current research in microelec- 
tronic devices is the development of a fast logic switch 
with minimal power dissipation per cycle. Devices based 
on spin wave interference 0] may provide an interesting 
alternative to conventional semiconductor gates by mini- 
mizing the need for current pulses. Recently, a spin wave 
NOT gate was demonstrated experimentally jlj . The de- 
vice consisted of a current-controlled phase shifter made 
by a ferromagnetic (FM) film on top of a copper wire. 
The application of a current along the wire creates a lo- 
cal magnetic field on the film, leading to a phase shift of 
its spin waves. 

In this letter we predict an effect that allows the de- 
sign of similar spin wave devices without the need for 
external current pulses or applied time-dependent mag- 
netic fields. We show that the dispersion of the lowest 
frequency spin-wave branch of a canted antiferromagnet 
depends strongly on the direction of spin wave propaga- 
tion. This occurs because of the long-ranged (dipolar) 
interactions of the magnetic excitations, which creates a 
gap for spin waves propagating with non-zero projection 
along the Neel axis. This effect allows electrical control 
of spin waves in multiferroic materials that possess simul- 
taneous ferroelectric (FE) and canted antiferromagnetic 
(AFM) order. 

Our model is applicable to the prominent room tem- 
perature muhiferroic BiFeOa (BFO) [1]. BFO films have 
homogeneous AFM order 0, Q , in contrast to the inho- 
mogeneous (cycloidal) AFM order present in bulk BFO 
0. The canted AFM order in BFO films is constrained 
to be in the plane perpendicular to the FE polarization 
P. Recently, Zhao et al. 0] demonstrated room tempera- 
ture switching of the Neel moment L = Mi — M2 in BFO 
films after the orientation of the ferroelectric moment was 
changed electrically. As we show here, spin wave prop- 
agation along P has high group velocity (^ 10^ cm/s), 
in contrast to spin wave propagation along L which has 
zero group velocity at fc = 0. Hence switching P for 
a fixed spin wave propagation direction allows electrical 
control of the spin wave dispersion, which assuming some 
loss rate will effectively stop long-wavelength spin waves 
such as those created in Q. 



Although a theory of AFM resonance for canted mag- 
nets was developed some time ago[l,Q, we are not aware 
of calculations of spin wave dispersion including magne- 
tostatic effects. The electromagnon spectra for a ferro- 
magnet with quadratic magnetoelectric coupling was dis- 
cussed without magnetostatic effects in Ref . [lO| and with 
magnetostatic effects in Ref. JT,. Recently we developed 
a theory of spin wave dispersion in bulk BFO, a cycloidal 
(inhomogeneous) multiferroic [l2l|. The lowest frequency 
spin wave mode was shown to depend sensitively on the 
P orientation because of the inhomogeneous nature of 
the antiferromagnetic order. Interestingly, we show here 
that BFO films with a homogeneous order display a sim- 
ilar effect, albeit due to a completely different physical 
reason: the magnetostatic effect. 

Our calculation is based on a dynamical Ginzburg- 
Landau theory for the coupled magnetic and ferroelectric 
orders. We assume a model free energy given by 
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+ {Jo + T]P'') Ml ■ M2 + dP ■ Ml X M2. (1) 

Here Mj is the magnetization of one of the two sub- 
lattices j = 1,2, and P is a ferroelectric polarization. 
The coordinate system is such that z points along one of 
the cubic (111) directions in BFO. The exchange inter- 
action J — ( Jo + rjP^) is assumed to have a quadratic 
dependence on P due to magnetostriction. The last con- 
tribution to Eq. Il]) is a Dzyaloshinskii-Moriya (DM) in- 
teraction, with a DM vector given by dP. Note that this 
changes sign under inversion symmetry, hence Eq. ^ is 
invariant under spatial inversion at a point in between 
the two sublattices. 

The design of multiferroic materials with enhanced 
couplings of this type was recently discussed [l3|- Al- 
though BiFeOa has no inversion center, its crystal struc- 
ture is quite close to an inversion-symmetric one, and the 
above free energy is derived by assuming that both the 
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FIG. 1: Spin and polarization waves in a canted multiferroic, 
such as a BiFeOa film. The sublattice magnetizations Mi, 
M2 lie in the plane perpendicular to the FE polarization P. 
Fluctuations SP denote polar phonons associated to vibra- 
tions of the FE moment, (a) Depicts the low frequency (soft) 
spin wave mode, (b) High-frequency (gapped) mode. The 
dots in the circle denote the position of the spins one quarter 
cycle later. The soft mode leaves the canting angle /3 invari- 
ant, while the gapped mode modulates /3. (c) Coordinate 
system. 



DM vector and polarization P are associated vifith the 
same distortion of the lattice. An alternative model for 
BiFeOa assumes the DM vector to be independent of P 
[3], i.e., requires Eq. ([T]) be invariant under spatial inver- 
sion at a point on top of one of the magnetic ions. Later 
we will discuss the implications of this alternative as- 
sumption for the electromagnon spectra, and show how 
optical experiments may determine which model is ap- 
propriate. 

The free energy is minimized by a homogeneous fer- 
roelectric and antiferromagnetic state, with FE moment 
(at E ^ 0) given by P = PqZ, with P^ = ^ + 0{(P). 
The magnetic moments are perpendicular to P, 

Mqi = Mo (sin (3x + cos (3y) , (2a) 
M02 = -Mo {- sin /3x + cos f3y), (2b) 

with canting angle (3 and magnetization Mq determined 
by tan/3 = (d-Po)/(J + J), and = ( J - r)/G, with 
j2 = (dPo)^ + J^- Below the Curie and Neel tempera- 
tures we have a < 0, and J > — r > respectively. 

Small oscillations away from the ground state are de- 
scribed by the Landau-Lifshitz equations. 
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where 7 is a gyromagnetic ratio. A corresponding set of 
equations is written for P in order describe the high fre- 
quency optical phonon spectra. Keeping only the lowest 
order in the deviations 5Mi and 5P, and focusing on the 
low requency magnetic oscillations we seek plane wave 



solutions of the type 

M, = Mq, + 5M,S''"-'^^\ P = Poz + SPe'^''-'-'^'\ 

(4) 

^From Eq. ([3]) we see that 6Mi must be perpendicular 
to Mi. Hence we may reduce the number of variables by 
using a parametrization for SMi shown in Fig. 1(c), with 
further definitions Y = yi + y2, Z = zi + Z2, y = yi — 1/2, 

Z = Zi- Z2. 

From Maxwell's equations we see that any macro- 
scopic wave producing nonzero fluctuations of SM ~ 
6M\ + bM2 must induce an AC magnetic h field. In 
the magnetostatic approximation this is obtained from 
V ■ h — — 47rV • 5M and V x /i « 0. The latter assumes 
the time variations are negligible in Maxwell's equations, 
which is a good approximation for spin waves provided 
k ^ wafm/c, with c the speed of light. For a canted 
AFM this is a good approximation provided the domain 
sizes are smaller than a few centimeters. The self induced 
field is therefore 



h = -Att {SM ■ fi) n. 



(5) 



where n is a propagation direction for the spin waves, 
fe — kh. The self-induced field contribute a term 
2tt{SM ■ nf to the free energy, tending to increase 
the spin wave frequencies whenever the quantity 5M = 
(— cos (/3)j/, sin {P)Y, Z) has a finite projection along h. 

In the magnetostatic approximation the linearized 
equations of motion are obtained by substituting 
Eqs. (21)- (IS]) into Eq. ([3]), and using the exphcit expres- 
sions for tan/3 and Mq. After some algebra the Landau- 
Lifshitz equations become 



-iuY + {J + J + ak^)Z -2h^ = 

ak'^Y + iujZ -2 sin f3hy = 

iCjz + {2 J + ak'^)y + 2 cos l3hx = 

{J — J + ak'^)z — iujy — 



-2d' cos f3S Pa:, (6a) 
-Ari' sin2f3SP,,{6h) 
-2d' cos2(35P^, (6c) 
-2d' sin (35Py, (6d) 



where we defined Co = lo/{'^Mq), d' = dMo, and 77' — 

Consider the pure spin waves in the limit SP —> 0. 
This case may be solved analytically, because the system 
of four equations decouples into two independent sets of 
equations on the variables {Y, Z) and (y, z). The former 
is a low frequency mode, because it corresponds to spin 
vibrations that leaves the canting angle (3 unchanged [the 
spins vibrate in phase, see Fig 1(a)]. The latter corre- 
sponds to spin vibrations half-cycle out of phase, leading 
to modulations of /3, and a high frequency gap equal to 
the DM interaction dPo [Fig- 1(b)]. Neglecting terms to 
second order in {dPo)/J, we may get an analytical ex- 
pression for the low frequency mode. 
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FIG. 2: (a) Low frequency magnetostatic spin wave dispersion 
for a BiFeOa film, for propagation angles 9 = Q° (propagation 
along the electric polarization direction i), 10°, 30°, 60°, 90° 
(propagation along the Neel direction y). The high frequency 
mode (not shown) has a gap equal to the Dyzyaloshinskii- 
Moriya coupling (5 x 10^° rad/s), and is nearly isotropic with 
respect to the direction of spin wave propagation, (b) Disper- 
sion including electrodynamical effects in the k < u)/c region. 
Note the relationship between the magnetostatic gap in (a) 
and the photon- magnon anticrossing in (b). 



This dispersion is anisotropic with respect to the po- 
larization (i) axis: For k along the x — z plane, we 
have a truly gapless mode to all orders in dPo/J, with 
u! « v2Jak. For k along y we find a gap equal to a frac- 
tion of the DM interaction, « ^/AitJJ {dP^) . This gap is 
a result of the magnetostatic correction in the presence 
of DM weak ferromagnetism. 

The physical origin of the magnetostatic gap is found 
by noting that 5M for a pure soft mode Y, Z Q,y, z — 
as fc ^ is approximately given by a rigid rotation 
around the z axis. In this limit, SM points exclusively 
along y, hence only propagation with some projection in 
this direction leads to a gap. 

A small anisotropy is also found for the high frequency 
mode {y,Zy^O,Y = Z = 0). For example, when fc || i; 
the high frequency mode gap increases to dPo y^l + iir/J. 

We calculated the coupled spin and polarization wave 
spectra solving the full set of Eqs. (pa|) - ((6d|) numerically, 
with parameters extracted from experiment Bi,!!- The 
low frequency spin wave branch within the magnetostatic 
approximation is shown in Fig. 2(a). The inset [Fig. 2(b)] 
shows the low frequency spectra beyond the magneto- 
static approximation, including electrodynamical correc- 
tions (For numerical convenience the speed of light was 
rescaled to 10^ cm/s). Note the anticrossing of the spin 
wave modes with the photon dispersion uj — ck, and the 
orientation dependence of the photon gap. As expected, 
we see that the strict k —^ limit has no orientation 
dependence. We emphasize that the latter low k limit 
is only observable for domain sizes of one cm or larger. 
The magnetostatic propagation anisotropy discussed in 



this work arises precisely because the spin waves travel 
with finite k > lu/c. 

Finally, we discuss the selection rules for the excitation 
or detection of magnon modes using an AC electric field. 
From inspecting the right hand side of Eqs. ([Sa|) and (|6b[) 
we see that the low frequency magnon may be excited 
electrically by the application of an AC field in the x 
or z direction. The former has a strong response in the 
presence of linear magnetoelectric effect {d 0), while 
the latter has a weak response (cx sin /3) due to magne- 
tostriction. 

The high frequency magnon {x,y) has a dielectric re- 
sponse only in the presence of the linear magnetoelectric 
effect, as seen in Eqs. (|6c|) and (j6d|) . This mode responds 
to electric fields in the y — z plane, with the z direction re- 
sponse larger by a factor of cos 2/?/ sin/3 ^ 2J/dPo » 1. 
The presence or absense thereof of this electromagnon 
using an optical probe may be used to discern whether 
the DM vector is linear in P as proposed e.g. in [l^ or 
if it is independent of P as suggested in 14 1. 

In conclusion, we predicted a magnetostatic gap 
anisotropy for the propagation of spin waves in a canted 
antiferromagnet. This effect may allow the electrical 
switching of magnons in multiferroic materials such as 
BiFeOa films. 
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